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Centers of near-infrared luminescence 
in bismuth-doped TlCl and Csl crystals 

V. O. SokolovQ V. G. Plotnichenko, and E. M. Dianov 

Fiber Optics Research Center of the Russian Academy of Sciences 
38 Vavilov Street, Moscow 119333, Russia 

A comparative first-principles study of possible bismuth-related centers in TlCl and Csl crystals 
is performed and the results of computer modeling are compared with the experimental data of 
Refs. [T}{3]. The calculated spectral properties of the bismuth centers suggest that the IR lumi- 
nescence observed in TlCliBi [l] is most likely caused by Bi^ • • • V^^ centers (Bi^ ion in thallium 
site and a negatively charged chlorine vacancy in the nearest anion site). On the contrary, Bi^ 
substitutional ions and BiJ dimers are most likely responsible for the IR luminescence observed in 
CsLBi |2113J. 

PACS numbers: 22.70.-a, 42.70.Hj, 78.55.-m, 



I. INTRODUCTION 

New near-IR luminescence, quite different from much- 
studied Bi^+-related luminescence [T, '5^, was discovered 
in 2001 in bismuth-doped aluminosilicate glass [6^. Then 
optical amplification at 1.3 pm was demonstrated [7]. 
Ever since bismuth-doped glasses and optical fibers based 
on such glasses attract a considerable interest due to a 
broadband IR luminescence in the range of 1.0 - 1.7 pm 
employed successfully in fiber lasers and amplifiers (see 
e.g. the review [8]). 

Although the origin of the IR luminescence is still not 
clear, recently a belief has been strengthened that subva- 
lent bismuth centers are responsible for the luminescence 
[9^. In our opinion, monovalent bismuth centers are of a 
particular interest. 

Glasses as disordered systems are very complicated to 
study impurity centers structure. In this regard, crystals 
may be of interest as model hosts containing bismuth- 
related centers. In particular, crystalline halides of mono- 
valent metals are convenient hosts for monovalent bis- 
muth centers. These crystals have a simple structure 
(primitive, Pm3m, or face-centered, Fm3m, cubic lat- 
tice). Bismuth can easily form monovalent substitutional 
centers in such lattice. Similar subvalent thallium and 
lead centers in Fm3m crystals were studied extensively 
(e.g. thallium in KCl [10^ and lead in MF2, M = Ca, Sr, 
Ba [H]). By analogy, the models of bismuth-related cen- 
ters in oxide glasses for fiber optics were suggested [12 . 

Bismuth-related IR luminescence in cubic halide crys- 
tals was studied for the first time in BaF2 [13], then in 
Csl (Pm3m) [2, 3^ and recently in TlCl (Pm3m) [l]. Two 
models of the centers in CsLBi were suggested in [2l|3], 
namely, a monovalent bismuth substitutional center, Bi^, 
and a dimer center, Bi2", formed by two Bi^ substitu- 
tional centers in the nearest cation sites with an extra 
electron. In what follows we report the results of com- 
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puter modeling of bismuth-related centers in TlChBi and 
CsLBi crystals. 



II. MODELING OF BISMUTH-RELATED 
CENTERS IN TICLBI AND CsI:BI CRYSTALS 

To study the origin of the IR luminescence in TlChBi 
and CsLBi, we performed a computer simulation of the 
structure and absorption spectra of three bismuth-related 
centers possibly occurring in both crystals. Basing on the 
assumptions made in [2j [3] and on analogy with Tl^(l) 
centers in alkaline halide crystals (see e.g. [10 ), we stud- 
ied the Bi^ substitutional center as the main form of 
bismuth embedding in TlCl and Csl, the Bi2" dimer cen- 
ter, and Bi • • • Vci and Bi • • • Vi complexes formed by the 
bismuth substitutional center and anion vacancy in its 
first coordination shell. 

The modeling was performed in a supercell approach. 
3x3x3 TlCl or Csl supercells (54 atoms) was chosen to 
model Bi^ and Bi • • • Vci or Bi • • • Vi centers, and 3x3x4 
supercell (72 atoms) was used for Bi2". In the central re- 
gion of the supercell certain cations were substituted by 
bismuth atoms and an anion vacancy was formed by a 
removal of one chlorine or iodine atom. Charged centers 
were simulated changing the total number of electrons 
in the supercell. Equilibrium configurations of bismuth 
centers were found by a complete optimization of the su- 
percell parameters and atomic positions with the gradi- 
ent method. All such calculations were performed using 
Quantum-Espresso package [14 in the plane wave ba- 
sis in the generalized gradient approximation of density 
functional theory with ultra-soft pseudopotentials built 
with PBE functional [15]. The pseudopotential sources 
were taken from the pslibrary v. 0.3.0 pseudopotential 
library [16 . 

To test the approach, TlCl and Csl lattice parame- 
ters were calculated for the unit cell and supercells with 
both atomic positions and cell parameters completely op- 
timized. The results convergence was tested with respect 
to the plane wave cutoff energy and to the k points grid. 




Figure 1. Calculated electron localization functions of Bi^ substitutional centers in l ( a)J TlCl and (b) 



Csl, and difference of 

calculated electron localization functions of Bi^ center from those of perfect lattice: (a) TlCl:Bi, (b) CsLBi (in the (111) plane). 



The energy cutoff > 700 eV and the number of k points with a relative accuracy better than 2% with only F point 

> 27 in the irreducible part of the unit cell Brillouin of the supercell taken into account and better than 1% 

zone were found to be enough to converge the total en- using 8 k points in the supercell in the irreducible part 

ergy within 10~^ eV per atom and to reproduce the ex- of the supercell Brillouin zone. The total energy conver- 

perimental lattice parameters with a relative accuracy of gence was not worse than that in the case of the unit cell. 
< 1%. The geometry of each supercell was reproduced 































o 


LD 
00 
O 


CO 
CZJ 


o 

CZJ 


LD 
O 





19.4 



14.3 
13.6 

11.8 
9.9 



21.9 



19.3 



13.6 



8.3 
7.9 



■22.6 



18.9 



14.7 



12.5 



9.2 



7.5 



4.2 





















00 
LO 

o 


o 


CM 

o 
































up 


?? 













23.9 



20.8 



17.3 



14.0 



6.2 



4.4 



25.0 



14.9 



12.5 



8.3 



6.1 





i 




1 






1 ' 




o 




i 










i 




LD 
O 


CN 
LD 
O 







23.0 



19.3 
18.0 

16.0 



10.1 
9.0 



(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



Figure 2. Calculated level and tran sition schemes of bismuth-r elate d centers: 
substitutional center in CsLBi 



(f ) Bi° • 



[(c)]Bi+ dimer center in TlCl:Bi,[(d) 
• Vi center in CsLBi. Level energies are given in 10*^ cm" 



BiJ dimer center in CsLBi 
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Bi^ substitutional center in TlChBi, 
Bi+ 
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Bi^ 



••Vci center in TlChBi, 



transition wavelengths in jim. 



Configurations of bismuth-related centers obtained by 
this means were used to calculate the electron localization 
functions using the programs from Quantum-Espresso 
package, the electron density distribution and effective 
charges of atoms by Bader's method using bader v. 0.28 
code [17 , and the absorption spectra of the centers by 
Bethe-Salpeter equation method based on all-electron 
full-potential linearized augmented-plane wave approach 
[T8| . The absorption spectra calculations were performed 
using Elk code [19 in the local spin density approxima- 
tion with PW-CA functional [20l|2T]. Spin-orbit interac- 
tion essential for bismuth-containing systems was taken 
into account. Scissor correction was applied to transition 
energies calculation with the scissor value found using 
modified Becke- Johnson exchange-correlation potential 
known to yield accurate band gaps in many solids [22H24] . 
The non-overlapping muffin-tin (MT) spheres of maximal 
possible radii R^^ were used. Convergence of the results 
was tested with respect to plane-wave cutoff energy, to 
the angular momentum cutoff for the MT density and po- 
tential, and to the k points grid choice. The plane-wave 
cutoff, kmax, was determined by the R^J^ ■ kmax = 7 
relation with R^J^ being the smallest MT radius. The 
angular momentum cut-off was taken to be / = 10. The 



self-consistent calculations were performed on the 3x3x3 
grid of k points uniformly distributed in the irreducible 
part of the supercell Brillouin zone. Further increasing 
the cutoff and k points density did not lead to significant 
changes in the results. The total energy self-consistence 
tolerance was taken to be 10~^ eV per atom. More dense 
4 X 4 X 4 /c points grid was applied to calculate dipole ma- 
trix elements in optical spectra calculations. 



Simplified configurational coordinate diagrams of 
bismuth-related centers were calculated in a model re- 
stricted to the lowest excited states with a displacement 
of bismuth atom(s) along [111] axis for Bi^ and Bi • • • Vj 
centers and along [001] axis for Bi2" center. In spite of 
the fact that the model is inherently rough, it shows that 
in all the centers studied the Stokes shift corresponding 
to a transition from the first excited state to the ground 
one do not exceed the accuracy of the excited state en- 
ergy calculation. Hence it seems reasonable enough to 
estimate the IR luminescence wavelengths by taking this 
Stokes shift to be zero. 
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Figure 3. Calculated electron localization functions of BiJ dimer centers in (a) TlC l and (b) Csl, and difference of calculated 
electron localization functions of BiJ dimer center from those of perfect lattice: (c) TlCliBi, (d) CsLBi (in the (101) plane). 



III. RESULTS AND DISCUSSION 

A. Bi^ substitutional centers 

Calculation of Bi^ substitutional center shows that 
both in TlCl and in Csl the crystal lattice is distorted 



rather slightly around it: bismuth atom lies in the cation 
site, the nearest chlorine or iodine atoms are displaced 
by almost 0.1a towards the bismuth cite, a being the 
crystal lattice constant (a ~ 0.3834 nm in TlCl and 
a ~ 0.4565 nm in Csl), and the nearest thallium or ce- 
sium atoms are displaced apart from the bismuth cite. So 



Bi-Cl and Bi-Tl distances are 0.3120 and 0.3889 nm, re- 
spectively (0.3320 and 0.3834 nm in perfect TlCl crystal) 
and Bi— I and Bi— Cs distances are 0.3641 and 0.4624 nm, 
respectively (0.3955 and 0.4567 nm in perfect Csl crys- 
tal). 

Bader analysis of the electron density around Bi^ cen- 
ter in TlCl showed the atomic effective charges to be 
-fl.03|e|, — 0.70|e|, and +0.73 |e| in -bismuth atom, in 
each of the nearest chlorine atoms, and in each of the 
nearest thallium atoms, respectively. The same analysis 
in Csl yielded +0.75 |e|, -0.67 |e|, and +0.80 |e| atomic 
effective charges in bismuth atom, in each of the nearest 
iodine atoms, and in each of the nearest cesium atoms, 
respectively. The effective atomic charges in perfect TlCl 
and Csl crystal lattices calculated by the same approach 
were found to be — 0.73 |e| and +0.73 |e| in chlorine and 
thallium atoms, respectively, and —0.80 |e| and +0.80 |e| 
in iodine and cesium atoms, respectively. 

In TlCl the total charge localized in cation cite and 
the neighboring anion sites was changed from —3.65 |e| in 
perfect lattice to —3.17 |e| in the lattice with bismuth sub- 
stitutional center. In Csl this total charge was changed 
from —4.00 |e| in perfect lattice to — 3.27|e| in the lat- 
tice with bismuth substitutional center. This shows that 
the electron density was displaced into the space between 
bismuth atom and the nearest anions and may be consid- 
ered as slight covalent contribution to Bi— CI and Bi— I 
interaction. In Csl such displacement is more noticeable 
(Fig.[l]). 

The calculated energy levels of the Bi^ center in 



TlChBi and CsLBi are shown in Figs. |2(a)| and |2(b)| re- 
spectively, together with the corresponding transitions. 
It should be noticed that splitting of the Bi^ ion states 
in bismuth substitutional centers in TlChBi and CsLBi 
is not described by crystal field theory due to total cu- 
bic symmetry of bismuth ion environment, unlike to the 
case of Tl^(l) center [10 . In this case the splitting is 
caused by electron density redistribution with a cova- 
lent contribution formed. Absorption near 1.0, 0.8, 0.7, 
and ^ 0.5 pm is found in the Bi^ center in TlCl. The 
above-mentioned excited states evaluative calculation al- 
lows one to expect a luminescence in TlChBi near 1.0 pm 
excited in this absorption bands. In CsLBi the IR lumi- 
nescence in the 1.2-1.3 pm range may be expected to 
be excited in absorption bands near 1.2, 0.7, 0.5, and 
> 0.45 pm. Besides, in TlCliBi another one lumines- 
cence transition, with much a lower lifetime, may occur 
near 0.8 pm. 



BiJ dimer centers 



respectively. The nearest chlorine or iodine atoms are 
displaced towards the dimer, and the nearest thallium or 
cesium atoms are displaced apart from the dimer. As a 
result, the Bi— CI and Bi— Tl distances are 0.3060 and 
0.3960 nm, respectively (0.3320 and 0.3834 nm, respec- 
tively, in perfect TlCl crystal), and the Bi— I, and Bi— Cs 
distances become -0.3314 and 0.4743 nm, respectively, 
as compared to 0.3955 and 0.4567 nm, respectively, in 
perfect Csl lattice. 

Bader's method analysis of electron density around the 
center shows that the excess electron charge, —1 |e|, is lo- 
calized almost completely in the first coordination shell 
of two bismuth sites, mainly in bismuth atoms and par- 
tially in the nearest chlorine or iodine atoms (Fig. |3|. In 
TlCliBi the effective charge in each of the bismuth atoms 
is found to be +0.82 |e|, and the effective charge in each 
of the nearest chlorine and thallium atoms is — 0.73 |e| 
and +0.72 |e|, respectively. So, as compared to Bi^ cen- 
ter, positive charge in each of the bismuth atoms is de- 
creased by 0.21 |e|, negative charge in each of the neigh- 
boring chlorine atoms is increased by 0.03 |e|, and positive 
charge in each of the nearest thallium atoms is decreased 
by 0.01 |e|. Hence the excess total charge localized in two 
bismuth atoms and their nearest neighbors turns out to 
be ^ — 0.6|e|. In CsLBi the effective charge in each of 
the bismuth atoms is +0.48 |e|, and the effective charge in 
each of the nearest iodine and cesium atoms is —0.70 |e| 
and +0.79 |e|, respectively. As compared to Bi^ center, 
positive charge in each of the bismuth atoms is decreased 
by 0.27 |e|, negative charge in each of the neighboring io- 
dine atoms is increased by 0.03 |e|, and positive charge in 
each of the nearest cesium atoms is decreased by 0.01 |e|. 
So the excess total charge in two bismuth atoms and their 
environment is ~ — 0.7|e|. So the center is actually Bi2" 
dimer both in TlChBi and in CsLBi. However, again a 
certain part of the electron density is displaced into the 
space between bismuth atoms and the nearest anions, as 
in Bi^ centers. 

The calculated energy levels of Bi2" center in TlCl 
and Csl and the corresponding transitions are shown in 
Figs. 2(c) and 2(d), respectively. Basing on the evalua- 
tive calculations of the first excited state one might ex- 
pect the IR luminescence bands in TlChBi near 1.1 and 
1.3 pm excited in absorption near 1.1, 0.8, 0.7, 0.5 and 
~ 0.4 pm, and in CsLBi near 1.2 and 1.6 pm excited in 
absorption near 1.2, 0.7, 0.6, and < 0.5 pm. One more 
luminescence band near ^ 2.4 pm in TlChBi and in the 
> 2.2 pm range in CsLBi corresponding to a transition 
from the lowest excited state might be observable at a 
low temperature. 



The modeling shows that the Bi2" dimer centers can 
actually occur both in Csl and TlCl crystals. 

In Bi2" centers the bismuth atoms are found to be dis- 
placed from the adjacent cation sites by 0.12a and 0.18a 
towards each other, so that the Bi— Bi distance is re- 
duced to 0.2903 and 0.2961 nm in TlChBi and CsLBi, 



C. Bi substitutional — anion vacancy complex 
centers 



The modeling of the Bi • • • Vanion complexes shows that 
the complexes in TlCl and Csl differ significantly in elec- 
tronic and spectral properties. The Bi • • • Vci complex 




Figure 4. Calculated electron localization functions of 
the (100) plane) a nd d ifference of calculated electron 



(a)]Bi+ • • • Vci center in TlCl:Bi and[(b)]Bi° • • • V? center in CsLBi (in 



ocalization functions of 



(a) 



perfect crystal and (b) Bi • • • Vj center in CsLBi from those of perfect crystal (in the (111) plane). 



Bi^ • • • Vq^ center in TlChBi from those of 



in TlChBi turns out to be similar in certain respect to 
Tl (1) centers in alkali halide crystals. However, the 
Bi • • • Vi complex in Csl differs strikingly from both the 
Bi---Vci andTl°(l) centers. 

In the Bi • • • Vci complex in TlCl the lattice relaxation 



turns out to be significant. The bismuth atom is dis- 
placed by 0.30a from cation site towards the vacant chlo- 
rine site, the nearest chlorine and thallium atoms are dis- 
placed towards bismuth atom, and thallium atoms sur- 
rounding chlorine vacancy are displaced apart from the 



vacant site. This results in Bi— CI, Bi— Tl, and Tl — Vci 
distances being 0.3024, 0.3598, and 0.4009 nm, respec- 
tively, as compared to the distances 0.3320, 0.3834, and 
0.3320 nm in perfect TlCl crystal. The relaxation is ac- 
companied by the electron density shifted from bismuth 
atom into the chlorine vacancy region, so that the com- 
plex center may be thought of as a bound pair o f ions , 
"Bi^ plus negatively charged V^^ vacancy" (Fig. 4(a)). 
Thus this center is a Bi^ • • • V^^ complex. 

In Bi • • • Vi complex in Csl, on the one hand, the lattice 
relaxation is as significant as that in Bi^ • • • V^^ center. 
Bismuth atom is displaced by 0.18a from the cation site 
towards the vacant iodine site, the nearest iodine atoms 
are displaced towards bismuth atom, but the nearest ce- 
sium atoms are displaced apart from bismuth atom. This 
results in Bi— Cs and Bi— I distances decreased to 0.4220 
and 0.3509 nm, respectively, from 0.4567 and 0.3955 nm, 
respectively, in perfect Csl crystal. On the other hand, 
the electron density distribution turns out to differ strik- 
ingly from that in Bi^ • • • V^^ complex in TlCl crystal. 
As distinct from TlCl, the electron density in Csl is not 
redistributed between bismuth atom and anion vacancy. 
So the Bi • • • Vi complex may be roughly thought of as 
a bound pair: a neutral substitutional atom, Bi , and a 
neutral iodine vacancy, Vj (Fig. 4(b)). In other words, 
this center turns out to be a Bi^ • • • Vj complex. 

Bader's method analysis of the electron density distri- 
bution around the center confirms these conclusions. In 
TlCl the effective charge of bismuth atom is +0.79 |e|, ef- 
fective charge of each of the nearest chlorine atoms is 
— 0.70|e|, and effective charge of each of the thallium 
atoms surrounding the chlorine vacancy is +0.66 |e|. So 
the excess (in comparison with perfect TlCl lattice) posi- 
tive charge localized in bismuth substitutional atom and 
its nearest neighbors is about 0.85 |e|, and the negative 
charge no less than 0.4 |e| is localized in the chlorine va- 
cancy and in its nearest neighbors. In Csl the effective 
charge of bismuth atom is — 0.02|e|, effective charge of 
each of the nearest iodine atoms is —0.82 |e|, and effective 
charge of each of the cesium atoms surrounding the iodine 
vacancy is +0.79 |e|. So the bismuth substitutional atom 
turns out to be practically a neutral atom, the excess (in 
comparison with perfect Csl lattice) positive charge lo- 
calized in its nearest neighbors is about 0.08 |e| and the 
positive charge no less than 0.06 |e| is localized in the 
iodine vacancy and in its nearest neighbors. 

The calculated energy levels of Bi^ 



Vqj and 



Bi • • • Vj cent ers a nd the co rresponding transitions are 
shown in Figs . |2 (e)] and [2(f)| respectively. Spectral prop- 
erties of the Bi^ • • • Vqj and Bi^ • • • Vj complexes may be 
understood in a model similar to Tl (1) center theory 
[10 . In such a model the complexes are considered as 
Bi^ ion or Bi^ atom, respectively, in the axial crystal 
field formed by neighboring chlorine or iodine vacancy. 
Obviously, the crystal field of a negatively changed chlo- 
rine vacancy is stronger than that of a neutral iodine 
vacancy. 

Three lowest states of Bi^ ion are known to arise from 



^P atomic state split by strong spin-orbital interaction 
in bismuth. The ground state of Bi^ ion is ^Pq. The 
first excited state, ^P^, is split by an axial crystal field 
in two levels, 6100 and 8300 cm~^, and the second ex- 
cited state, ^P2, is split in three levels, 12500, 14900 and 
25000 cm~^. In a free Bi^ ion the electric dipole (El) 
transitions between three spin-orbital components of the 
^P state are forbidden but under the infiuence of the 
crystal field the transitions become allowed. Basing on 
the above-mentioned evaluative calculations one expects 
the IR luminescence in two bands near 1.6 and 1.2 pm, 
both excited in absorption near 0.8, 0.7, and ~ 0.4 pm. 
A luminescence band with a significantly shorter (by an 
order of magnitude) lifetime may occur near 0.8 pm. 

The ground state of Bi atom is known to be ^S3/2. 
The first excited state, ^03^2, is split by an axial crys- 
tal field in two levels, 9000 and 10100 cm~^. The sec- 
ond excited state, ^05/2, is split in three levels, 16000, 
18000, and 19300 cm~^, and the forth excited state, 
^Pi/2 5 is not split by an electrostatic field. Notice that 
the splitting is expected to be considerably slighter than 
that in Bi^ • • • V^^ complex. El transitions from the 
ground state of Bi • • • Vj center, corresponding to the 
^^3/2 atomic state, to all the states arising from the ^D 
one, turn out to be weak since in a free atom such tran- 
sitions are parity-forbidden. Hence the only relatively 
intensive absorption band near 0.43 pm corresponding to 
^83/2 -^ ^Pi/2 transition may be expected to occur in the 



Bi^ • • • Vj center (Fig. 2(f)). With the above-mentioned 



evaluative calculations taken into account one might ex- 
pect that 1.0 - 1.1 pm luminescence corresponding to 



D3/2 

band. 



53/2 transitions is excited in this absorption 



IV. CONCLUSIONS 

Spectroscopic data [THS] and the present results of 
modeling of bismuth-related centers in TlChBi and 
CsLBi crystals suggest that different centers are respon- 
sible for the near-infrared luminescence in these crystals. 

In TlCliBi the luminescence observed near 1.2 pm is 
caused mainly by Bi • • • Vci complexes formed by Bi^ 
substitutional ions and negatively charged chlorine va- 
cancies. However, the bismuth-related contribution to 
the total absorption is caused mainly by single Bi^ sub- 
stitutional centers not contributing to the IR lumines- 
cence. This explains a significant difference between the 
IR luminescence excitation spectrum and TlChBi absorp- 
tion spectra [1 . Bi^ substitutional centers may cause a 
luminescence near 1.0 pm not observed in y[j. 

In CsLBi the IR luminescence is caused mainly by two 
types of bismuth centers: Bi^ substitutional centers give 
rise to a luminescence observed in the 1.2 - 1.3 pm range, 
and Bi2" dimer centers cause a luminescence observed 
near 1.2 pm and near 1.6 pm. These conclusions agree 
with the experimental data of [2l, ^ and confirm the as- 



sumptions made there. Besides, dimer centers may cause 
possible luminescence in the ^2.2 pm range not observed 
inRefs. 0[3]. 

Bi2" dimer complexes in TlChBi can contribute percep- 
tibly neither to the IR luminescence spectra nor to the 
absorption spectra, as distinct from CsLBi crystals. 

Bi • • • Vi complexes formed in CsLBi by Bi substitu- 
tional atoms and neutral iodine vacancies may contribute 
to the IR luminescence only near 1.0 pm, as distinct from 
TlCliBi crystals. One might suppose that the lumines- 
cence in the 0.95-1.15 pm range observed in [2 at a low 
temperature under 0.4-0.5 pm excitation is contributed 



by Bi^ • • • Vj centers formed due to electrons capturing 
in Bi^ substitutional centers in the vicinity of iodine va- 
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